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Abstract. We show that Zipf's Law for the largest clusters is not valid in an 
exact sense at the critical point of the fragmentation phase transition, contrary 
to previous claims. Instead, the extracted distributions of the largest clusters 
reflects the choice of universality class through the value of the critical exponent 
r. 
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1. Multifragmentation 

Nuclear multifragmentation offers a unique opportunity to study phase transitions 
under extreme finite size constraints. Throughout the last decade models of this 
phase transition have been developed and tested against experimental data. In 
particular, recent scaling analysis of multifragmcnt emission has yielded important 
information on the character of the phase transition, yielding values for the critical 
temperature and the values of the exponents r and a E]- 

Event samples provided by multifragmentation are also a rich source for sta- 
tistical analysis. As one example, we cite the observation that intermediate-mass 
fragment multiplicity distributions as a function of the total transverse energy fol- 
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low sub-Poissonian distributions [ 3 . Our group, however, was able to show that 
this binomial scaling was exclusively a finite-size effect El El- Furthermore, we 
were able to show that the same patterns exist in word-length frequencies in the 
literature, where the system size is given by the length of a paragraph [[7]. 

In 1949, the linguist G. Zipf published [ 8 another very interesting observation 
about many languages: if one ranks all words in a wide variety of texts in a given 
language, then the rank is inversely proportional to the frequency of occurrence 

F{r) = cr~^ (1) 

(with c = constant. Awl). The same distribution was found much earlier by Pareto 
in ranking companies by their income [ El ■ Similarly, ranking US cities by their 
population yields the same law. Recently, Watanabe [ TU observed in numerical 
simulations that percolation clusters at the critical threshold also obey Zipf 's Law, 
and Ma [ 1111 IT^ proposed to make use of this finding to detect a crossing of the 
critical threshold. This motivates the current study. 

2. Cluster Sizes at the Critical Point 

In a scaling theory the cluster size distribution follows the scaling function 

n{A,e) = aA-^f{A- -e) (2) 

Here e is the fractional deviation of the control parameter from the critical value, 
e = {p — Pc)/Pc, A is the cluster size (— number of nucleons in the case of nuclear 
fragmentation), / is the scaling function with the boundary condition of /(O) = 1, 
and a and t are the critical exponents that determine the universality class of the 
phase transition. The normalization constant a can be obtained from the condition 
that all nucleons have to belong to some cluster: 

V 

Y,A-7i{A,e)^V, (3) 

A=l 

where V is the event size {— total number of nucleons). 
Let us first address scaling at the critical point, 

e = ^ n{A, 0) = N{A) = a A-\ (4) 

i.e. the cluster size distribution is a pure power law. 

The probability Pist{A) that a cluster of size A is the largest in a given frag- 
mentation event is the product of the probability that there is at least 1 cluster of 
size A present, and that there are clusters of size larger than A present in the 
same event. 



Pist{A)=p<i{A)-p^{>A) = [l-po{A)]-po{>A). 



(5) 
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Since the average yield of clusters of size A is N{A) = aA , we obtain for the 
average yield of clusters of size bigger than A the sum 

V V 

N{>A)^ N{t)= at-- = aC{T,l + A)~aar,l + V) (6) 

i=A+l i=A+l 

where C is the generalized Riemann function, and we find for the normalization 
constant 

V 

a^V/YA'-^ -V/Hl-- (7) 

with Hy-"^ is the V*''^ harmonic number of order 1 — t. 

For V ^ A, the probabilities to have k clusters of size ^ in a given event are 
Poissonian,'^ 

Pk{t) = l(7V(z))'=e<^«> (8) 

Thus, we finally obtain for the probability that a cluster of size A is the largest 
in a given fragmentation event 

P^sM) = [1 - Po{A)] ■ po{>A) = [1 ~ e--""'^] ■ e[<(-^i+^)-<Ki+^)l, (9) 

with a given by equation |7| 

From this we obtain the exact expression for the average size of the largest 
cluster by summation over all A: 

V 

(Ai,t) = ^^-Pi,t(A) (10) 

The probability that a cluster of size A is the second biggest in a given event 
is the sum of two terms. The first term is the probability that two or more clusters 
of size A are present times the probability that none bigger than A is present. And 
the second term is the probability that exactly one cluster of size larger than A is 
present times the probability that at least one cluster of size A is present. This 
yields: 

P2nd{A) = p>2{A) ■ Pq{> A) + p>i{A) ■ pi{> A) 

= [l-po{A)-p,{A)]-po{>A) + [l-po{A)]-p,{>A) (11) 

Obviously, this can be continued via recurrence relations to the third, fourth, ... 
nth biggest clusters. It is important to note that these exact analytical results only 
depend on only two parameters, the size of the fragmenting system, V, and the 
value of the critical exponent r. 

As an example, in figure ^ we show the probability distributions for the largest 
through fourth-largest fragments as a function of mass number A for two different 
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Fig. 1. Probability distributions for the largest to fourth- largest clusters in a 
fragmentation event of a system of size V = 100 (top) and V = 10,000 (bottom). 
The arrows indicate the corresponding average values. 



system sizes, V ~ 100 and V = 10, 000. Also marked in this figure, by the vertical 
arrows, are the average values for the largest through fourth-largest fragment. In 
both cases the value for r was chosen as 2.18, that of a 3d percolation universality 
class. Next we investigate whether Zipf's Law, equation ^ is satisfied. 
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Fig. 2. Average size of the largest fragment of rank r for a system size of ^ = 
10, 000. The solid lines are power-law fits. The different sets of data were generated 
with different values of the critical exponent r, as indicated in the inset. 

In figure 121 we plot the average size of clusters of rank r for 5 different values of 
the critical exponent r in a system of size V = 10, 000 (plot symbols), together with 
their best respective power law fits of the form {Ar) oc r^'*'. It is clear that these 
data approximately follow power laws, but not exactly. The power law exponents A 
numerically extracted are 1.26, 1.16, 1.07, 0.97, 0.7 for r = 2, 2.18, 2.3, 2.5, 3, respec- 
tively. One can see that in all cases our results are close to Zipf's Law. But there 
is a systematic dependence of A on the critical exponent r, which is something that 
should have been expected. 

If one uses an expansion for our analytical results, one finds that the average 
cluster size as a function of rank r follows a more general Zipf-Mandelbrot distri- 
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bution [nni ini 

{Arth) = c[r + k)-^ (12) 

where the offset k is an additional constant that one has to introduce, and A is 
asymptoticahy approximated as a function of the critical exponent t 

A(r) = (13) 

T — 1 

One obtains almost perfect fits to the average cluster sizes as function of their 
rank by using equation El This is particularly evident if one plots the ratio of the 
average size of the largest cluster to that of the rth largest cluster, see left side of 
figureOl The resulting value of the exponent A as a function of the critical exponent 
T are shown on the right side of figure 13 The dashed line is the approximation of 
equation 1131 
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Fig. 3. Left: Ratio of the average size of the largest cluster to that of the rth 
largest. Dots are the analytical results, thin lines the best simple power law fits, 
thick lines the fits with the distribution of equation 1121 Right: Effective exponent 
A of the Zipf-Mandelbrot distribution fit to the exact fragmentation cluster size 
distribution as a function of the critical exponent r for three different system size, 
V = 100, V = 10,000, and V = 1,000,000. Dashed line: asymptotic expansion of 
equation 1131 



3. Conclusions 

We have shown that the distribution of largest clusters does not exactly follow the 
functional form implied by Zipf's Law, but instead a more general distribution of 
the Zipf-Mandelbrot kind, and with a power-law exponent that is not unity, but 
varies as a function of the critical exponent t. Thus the observation of frequency 
distributions close to the expectation of Zipf's Law does not yield additional physics 
insight and must be considered coincidental. 
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d. We have also conducted studies in the percolation model, where the independent 
Poissonian emission condition docs not strictly hold, with similar results. 

e. After this project was completed, we became aware through private discussion 
that X. Campi and H. Krivine have arrived at somewhat similar conclusions 
independently, and in a different way [I15|. 
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